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(54) Title: EPILEPTIC SEIZURE DETECTION BY NONLINEAR MEraODS 
(57) Abstract 

Methods and apparatus for automatically detecting epileptic 
seizures by monitoring and analyzing brain wave (EEG or MEG) signals 
of a patient (11) using EEG electiodes (13). Steps inchide: acquiring 
the brain wave data from the patient; digitizing the data; obtaining 
nonlinear measures (15) of the data via chaotic time series analystr, 
obtaining lime serial trends in the nonlinear measures; determining that 
one or more trends in the nonlinear measures indicate a seizure, and 
providing notification of seizure ocmrrence. 
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EPILEPTIC SEIZURE DETECTION BY NON-LINEAR METHODS 



The United States Government has rights in this invention puisuant to contract no. DE- 

AC05-84OR21400 bet>wen the United States Departmem of Energy and Loctteed 
Energy Systems, Inc. 

5 FIELD OF THE fNVENTION 

The present invention relates to the application of chaotic time series analysis (GTS A) to 
electroencephalogram (EEG) data and magnetoencephalogram (MEG) data, and more 
particularly to the analysis of the data to detect qrileptic seizures. 

CROSS-REFERENCE TO RELATED APPUCATIONS 
10 niis invention is related to Docket No. I743-X EpUeptic Seizure Prediction by Non- 

linear Methods by Lee M. Hively. Ned E. Clapp. C. Stuart Daw. and WiUiam F. Lawkins and to 
Docket No. 1 744-X Method and Apparatus for Extraction of Low-Frequency Artifects from 
Brain Waves for Alertness Detection by Ned E. Clapp and Lee M. Hively, both of which are 
filed on even date herewith, and both of which are assigned to the same entity. 

15 BACKGROUND OF THE INVENTION 

TTie theory of nonlinear dynamics provides a basis for understanding and potentialjy 
controlling many complex physical and engineering systems. An extensive literature exists for 
nonlinear dynamics in the brain and related work (1 8). It is weU known that biain waves exhibit 
seemingly random, unpredictable behavior, that is characterisdc of detenninistic chaos (1 , 20). 
20 Moreover, chaotic behavior is "normal,- whUe nonchaotic or periodic behavior b indicative of 
pathophysiology in experimental epUepsy (20). Schifif et al. (1 1) showed that chcmicaUy- 
induced seizures in lat-brain can be electrically controUed. leading to specubtion (4) that human 
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epilq3sy may be controlled without drug or surgical intervention. However, effective use of 
chaos control for epilepsy requires definitive seizure detection. Thtis, this invention diagnoses 
brain wave data via chaotic time series analysis (CTSA) methods to detect an epileptic seizure. 

Nonlinear analysis of neurological diseases via EEG data is extensive. For example, see 
the 1994 review by Elbert et al. (18). Epilepsy can be recognized only with clear EEG 
manifestations, but even these seizures are not easy to detect because there is no stereotyped 
pattern characteristic of all seizures (S). Woric by Olsen and colleagues (7) used various linear 
measures with autoregressive modeling, discriminant analysis, clustering, and artificial neural 
networks. Valuable nonlinear tools for studying EEG data include correlation dimension, mutual 
information fimction, Kolmogorov entropy, phase-space attractors, and largest Lyapunov 
exponenL 

Very recent analysis by Theiler (16) stiutied correlation dimrasion and Lyapunov 
exponent, using a form of surrogate analysis on asingle EEG time series during an epileptic 
seizure. The surrogate analysis involved a raxidom shufiling of blocks of time serial data, each 
block containing one quasi-periodic spike*wave complex. The auto-correladon fimcdon for the 
original data is nearly indistinguishable bom the surrogate data. The correlation dimension for 
the original data is significantly different from the surrogate data only at large scale sizes and 
large embedding dimensions. The maximum Lyapunov exponent (X) was negative for both the 
original and surrogate data and not substantially different, contrary to previous work which found 
positive k values. Theiler concluded that his analysis suggests a nonlinear oscillator with noise 
on the time scale of the spike-wave complex, but cannot indicate whether chaos exists on a 
shorter time scale. 

Other patents for epileptic seizure detection have been granted. U. S. Patent No. 
531 1,876 "Automatic Detection of Seizures Using Electroencephalogr^ic Signals'" by D. E. 
Olson et al. relies on linear analysis of EEG features for seizure detection. U. S. Patent No. 
5,349,962 "Method and Apparatus for Detecting Epileptic Seizures" by J. S. Lockard et aL 
detects seizures if Aewavefom is withm a predeterinined threshold. These patents use linear 
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methods, as distinguished from the nonlinear techniques of the present invention. No other work 
is known that applies several CTSA measures to EEG data or MEG data for systqnatic 
chaiactoization of non-seizure, seiane. and transition-to.sei2uie. 
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OBJECTS OF THE INVENTION 
Accordingly, it is an object of the present invention to provide new and improved 
methods and apparattis for detecting epileptic seizures in a patient and providing notification to 
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pennit assistance to be given to the patient or to a person who can assist the patient 

Further and other objects of the present invention wUl become apparent from the 
description contained herein. 

SUMMARY OF THE INVENTION 

In accoidance with one aspect of the present invention, the foregoing and other objects 
are achieved by a method for automaticaUy detecting an epileptic seizure in a patient comprising 
the steps of: providing at least one chamiei of a patienf s raw brain wave data, called e-data, 
selected from the groiq, consisting of electroencephalogram data and magnetoencephalogram 
data; separating the e^ata into artifact data, caUed f^ and artifict-fiee data, called gHlata, 
while preventing phase distortions in the data; processings^ through a low-pass filter to 
producealow-pass-ffltered version of g.data.caUedhKlata; applying at least one noUinear 
measure selected from the group consisting of time steps per cycle. Kohnogorov entropy, first 
minimum in the muaial information function, and correlation dimension to at least one type of 
data selected from the group consisting of e-data, f-data. g^ and h^laia to provide at least 
time serial sequence of nonlinear measures, from which at least one indicative trend selected 
fom the group consisting of abnipt increases and abnipt decreases can be deters 
comparing at least one indicative trend with at least one known seizure indicator, and 
detemuning from said comparison whether an cpUeptic seizure is occurring in the patient 

In accordance with a second aspect of the present invention, the foregoing and other 
objects are achieved by apparatus for automatically detecting an epUeptic seizure in a patient 
comprising: data provision means for providing at least one chamiel of a patient's raw brain 
wave data called e-data selected from the groiq, consisting of electroencephalogram data and 
magnctoencephalogramdata; separation means for separating e-data into artifact data. caUedf- 
data. and artifect-free data. caUed g-data, while preventmg phase distortions in the data, 
communicably comiected to the data provision means; low-pass filter means for filtering g-data 
to produce a low-pass filtered version of g-data, called hslata, communicably comiected to the 
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sq>aration means; application means for applying at least one nonlinear measure selected from 
the group consisting of time steps per cycle, Kolmogorov entropy, first minimum in the mutual 
infoimation function, and conelation dimension to at least one type of data selected from the 
group consisting of enlata, f-data, g*data, and b-data to provide at least one time serial sequence 
of nonlinear measures, from v^cb at least one indicative trend selected from the group 
consisting of abrupt increases and abrtipt decreases can be determined, communicably connected 
to the low*pass filter means; comparison means for comparing at least one indicative trend with 
known siezure indicators, communicably connected to the application means; and determination 
means for determining from the comparison whether an epileptic seizure is occurring in the 
patient, communicably connected to the comparison means. 

BRIEF DESCRIPTION OF THE DRAWINGS 
In the drawings: 

Fig. 1 is a block diagram showing how, in accordance with the present invention, EEC 
data is obtained from the patient, digitized, processed, and analyzed by nonlinear methods to 
detect and/or predict epileptic seizures. 

Fig. 2 shows standard EEG electrode positions on a patient's scalp for the bipolar 
montage, looking from above. 

Fig. 3 shows sample plots of EEG data, to illustrate raw data (c-data as the dotted curve 
in fig. 3a), and artifect-filtered data (g-data, shown in fig. 3b), as an example of the method and 
apparatus of this invention. 

Fig. 4 comprises Figs. 4.1, 42, 4.3, and 4.4, and shows linear and nonlinear measures of 
time serial data for Example I. 4.1 shows linear and nonlinear measures of raw (c-data); 4.2 
shows linear and nonlinear measures of artifact (f-data); 4.3 shows linear and nonlinear measures 
of artifect-filtered (g-data); and 4.4 shows linear and nonlinear measures of artifact- and low-pass 
filtered (h-data). Various ma&ematical properties and characteristics arc shown for each data 
type, and are further shown in parts (a), (b), (c), (d). (e), (f), (g), (h), and 0) for each data type. 
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Tlicse measures are shown as curves, each curve representing a time serial sequence of linear or 

nonlinear measures. Within each curve, significant features such as abrupt increases and abrupt 

decreases may be viewed as indicative trends which are then compared to to trends which have 

been shown to be siezure indicators. Thus from the comparison it can be detennined whether a 

5 siezure is occurring in the patient Inpaits (a), (b).(cX(d). and (e) for each data type, the soUd 

line is the specific measure, the dashed line is the 1 I^wim average, and the dotted line is tte 
standard deviation for the measure. 

Fig. 5 comprises Figs. 5. 1. 52. S3, and 5.4. and shows linear and nonlinear measures of 
time serial data for Example H. 5.1 shows linear and nonlinear measures of raw (e-data); 52 

10 shows linear and nonlinear measures of artifiwt EEC (f-data); 5.3 shows linear and nonlinear 
measures of arti&ct-filtered (g^); and 5.4 shows linear and nonlinear measures of artifact- 
and low^ filtered Owlata). Various matiiematical properties and characteristics as computed 
for each data type are fiirtiier shown in parts (a), (b), (c). (d). (e), (f). (g), (h). and © for each data 
type. These measures are sho>yn as curves, each curve representing a time serial sequence of 

15 linear or nonlinear measures. Within each curve, significant featims such as abrupt increases 
and abrupt decreases may be viewed as indicative trends which are then compared to to trends 
which have been shown to be siezure indicators, nius from the comparison it can be determined 
whetiier a siezure is occurring in tiie patient In parts (a), (b), (c), (d), and (e) for each data type, 
the soUd line is die specific measure, the dashed line is the 1 1 -point average, and die dotted line 

20 is the standard deviation for the measure. 

For a better uaderstaoding of the present invention, togetiier with odier and finther 
objects, advantages and capabilities thereof, reference is made to the foUowing disclosure and 
appended claims in connection with the above-desciibed drawings. 



'25 



data. 



DETAILED DESCRIPTION OF THE INVENTION 
Chaotic time series analysis (CTSA) is applied to human electroencephalogram (EEG) 
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For reference, equation manbers are shown at the right of each equation. 

TlBee epochs were examined: epileptic seizure, non-seizure, and transition 6om non- 
seizure to seizure. The CTSA tools were applied to four forms of these data: raw EEG data (c- 
data), artifect data (f-data) via application of a quadratic zero-phase filter of flic raw data, 
artifect-filtcrcd data (g-data) that was the residual after subtracting f-data firom e-data, and a low- 
pass-filtered version (h-data) of g-data. Several nonlinear measures umquely indicate an 
epileptic seizure, including: an abrupt increase in the number of timesteps per cycle for f-data; 
an abrupt increase in the number of timesteps per cycle for g-data; an abrupt increase in the 
number of timesteps per cycle for h-data; an abrupt increase in the entropy of e-data; an abrupt 
increase in the entropy of f-data; an abrupt increase in the entropy of g-data; an abrupt increase 
in the entropy of h-data; an abnipt decrease in the first minimum in the mutual infoimation 
function for f-data; an abrupt decrease m the first minimum in the mutual infonnation function 
forb-data; an abrupt increase in the correlation dimension of e-data; an abrupt increase in the 
correlation dimension of f-data; an abrupt increase in the correlation dimension of g-data; an 
abrupt increase in the correlation dimension of h-data; and combinations thereof! Analysis of e- 
data shows that statisticaUy significant nonlinear strucoire is presem during the non-seizure, 

tranation, and seizure epochs. 

Two sets of channel-l3 EEG data from one patient are provided as examples, to Fig. 1 of 
the drawings. 11 shows the patient's head, looking from above. 13 shows EEG electrode 
positions on the patients scalp. 15 shows a nonlinear measure of EEG data. In Fig. 2 of the 
drawings, C13 labels tiie position where tiie channel 13 data, whidi is used in tills woik, 
originates. Analysis of C13 data was chosen for these examples to demonstrate die robust 
removal of cye-blink artifact, which otiierwise dominates channel 13 because of its proximity to 
tiieeye. TTus mcdiod and apparams can be appUed to data from other EEG channels, as well as 
to MEG dala.asbapparenttotiK)se skilled in tiieart Botii sets of data included non-seizure, 
transition-to-seizmc, and epUeptic seizure data. TTie analysis included various linear measures 
(standard deviation, absolute average deviation, skewedness. kurtosis), plus nonlinear measures 
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(time steps per cycle, Kolmogorov entropy, fim mimmum in the mutual infotmation function, 
and coirelation dimension). Four fonns of the data were analyzed: raw EEG (e) data, aiti&ct (f) 
data via application of a zero-phase quadratic filter, arti&et.filtered (g) data that was the residual 
after subtracting fsiata fiom g-data, and a low-pass-filtered version (h) of the gnJata. TTie 
nonlinear measures clearly discerned the epUeptic seizures in these exanq>les. whUe none of the 
linear measures provided definitive seizure indication. Nonlinear tools identified seizure 
indicatora with and without aiti&ct removal, with and without low^ filtering, demonstrating 

the robustness ofthese methods to noise and artifiuas. Surrogate analysis of e-data showed that 
this data has significant nonlinear structure. 

In the examples described herein, sixteen chaonels of EEG data were analyzed in the 
bipolar montage, as illustrated in Fig. 2. Hie data were retrieved in analog form from VHS tapes 
and converted to digital form with 12-bit precision, giving an integer between -2048 and +2047. 
The digital sampling rate (f J was 5 1 2 Hz over a total sample time of 1 0-23 minutes, 
corresponding to a total dataset size of 9.8-22.5 megabytes in binary form, nuee epochs of data 
15 were examined: epUeptic seizure and post-seizure. non-seizure, and transition fitnn non-seizure 
state to seizure (transition). 

It is acknowledged that to detect an epileptic seizure as it occurs, the brain v^^^ 
would not be recorded data, but would be currcnUy-occurring data. TTiis data would be taken 
fiom the patient to the apparatus direcdy using standard EEG or MEG methods, or indirectly by 
transmitting the data to an apparatus remote fiom the patient by means such as telephone, radio, 
or other communications means well known to the skilled artisan. 

Raw brain wave data contains not only signals associated with brain activity, but also 
contains arti&cts (e.g., eye blinks, muscle twitches, chewing, etc.) that obscure the brain-wave 
signal. In order to observe artifiict-frec data and artifact data independenUy, the raw data must be 
separated into artifett data and artifect-free data. A zeio-phase filter was developed and used to 
remove low-frequency artifects. based on the foUowing criterion. A zero-phase-shift filter was 
needed to prevent phase distortions when subtracting the filter output (the "artifiict" data signal) 
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from the raw data signal to yield an undistorted artifact-filtered or artifact-free data signal, 
because phase relationships are most important in the subsequent nonlinear analysis. Standard 
high-pass filter techniques do not meet tiiis criterion. A computationally fast, simple, low- 
frequency signal follower was necessary to eventually apply the filter in real- or near-real time. 
Consequendy, quadratic regression analysis was used, with the same number of data samples on 
eidier side of a central point Other standard digital filtering methods (IS) could not meet this 
requiremenL 

The zero-phase filter method and apparatus v*ich may be embodied in various ways well 
known to one skilled in the art, such as a specially-programmed computer or a programmed 
integrated circuit, semi-conductort chip, or micro-processor, is as follows. For a ^)ccific 
channel, the EEG signal (e) at time (t) is sampled at regular intervals ft =i At) to yield a set of 
time serial data ei«e(tO. We choose a filter-window length of 2n+ 1 points from the time series, 
where n is the number of points on either side of the central point (ej as indicated in the 
sequence below. 

central point 
I 

n points n points 

The data was fitted to a quadratic equation that takes the form: Fj * ¥{Xd = a, ft - tj^ + aj 
tt-tJ + aj^aiTi' + ajTj + aa. Here, te = c At is the time at the central point, and = tj - 1^. This 
approximation is fitted to the data, by minimizing the sum of squares of the differences between 
the quadratic equation, F(t), and the raw EEG data, e(t), corresponding to the minimum in the 
following fimction: 
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The minimum in L is found from the condition dUda, - 0. for k - { 1 . 2. 3 } . fonning 
three simultaneous linear equations in three unknowns. TTie wndow^yeraged artifact (FJ is 
given by the fitted value of the centHd point. F. - F(0) » a. Note that the sums o 
ofT, are zero and tiat symmetric sums over even powers of T, (over i fiom -n to can be 
5 converted to sums from 1 to n with T, = i At. yielding a window^vcraged solution for ti« artifect 
signal: 

(4«^+4»-3)(2»+l) 

Here, Ei indicates tiie sum over i fiom * to +n. Sums over even powers of "i" were 
explicitiyevaluatedwithstandardformulae(lO). The efifort to evaluate F. can be reduced 

10 substantially by computing the sums initiaUyfiomEq^ (at c = n+l).and then usingUie 
following recursions thereafter 



S ' •■«c^*ir«,.r Eq.3 



■ 



n 

TTic light-haod sides of Eqs. 3-5 only involve the sums previously computed. 
AppUcation of Eqs. 2-5 to the N-point set of original time serial EEG data (e,. illustrated as the 
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dotted curve in Fig. 3a) yields an artifact dataset (fj or f-data, illustrated as the solid curve in Fig. 
3a) with (N-2n) points that contains the low frequency artifact signal. The residual signal (g^ or 
g-data, illustrated as the solid curve in Fig. 3b) is the difference, & Cj - , and is a signal that is 
free of low-frequency artifects. Subsequently, the g-data is processed through a standard fourth- 
order low-pass filter at 50 Hz (e.g., see Ref 9) to yield artifact-filtered, low-pass-filtered data (hi 
or h -data) that is free ofboth low- and high-frequency artifacts. Note that spike-wave 
phenomena at 100 Hz in h-data are attenuated by 28 db (a factor of 25), while the g-data retain 
the fiill spike-wave signals. A standard second-order, third-order, or fourth order filter at 
frequencies between about 35 Hz and about 60 Hz would also be effective as a low-pass filter. 

The filter-window lengdi (n=128) corresponds to a frequency of 2.0 Hz [=512 Hz/(2n 
+1)]. Figure 3a shows an example of the application of this method, with (raw) c-dala in light 
gray and a superimposed (dark line) artifect signal (f-data), which clearly follows the low- 
frequency trends. Figure 3b shows the residual signal (g^iata) for this example, as having Utf^ 
low-frequency component while retaining the higher frequency information. 

For a specific EEG channel, a time history of &e nonlinear measures was obtained by 
applying the CTSA tools to a scries of 20-second analysis-windows of the four data types (e, f; 
g. h). These data are designated herein as x^. The length of the analysis window (w) was 10^40 
points. Each analysis-window had a 5,l20.point overlap with the previous (or next) analysis- 
window of data. This 50% overlap provides an optimal mix of new and old data for smooth 
time-history trend generation (2), as illustrated below. 

|<--0 to 20 seconds— >| 

|<:~10 to 30 seconds— >| 

|<— 20 to 40 seconds— >| 

|<-100 to 120 seconds->| 

|<-1 10 to 130 seconds->| 
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10 



15 



20 



25 



TTie zero-phase quadratic filter provides arti£u:t-filteced data with fiequencies of i2 Hz. 
A heuristic for (linear) Fourier analysis is that ^ 10 periods of data are required to feithfully 
recover cycUc information at a specific frequency. Tlus. ^5 seconds ofdata are needed to obtain 
FourierampUtudeandphaseinfomuuionatasignal frequency of2H2L However, this heuristic 
doesnotapplytononlinearanalysis. For example. 20 seconds of data are necessary to obtain 
consistemresuItsfortheKohnogoiDvemropy. TTiis need for longer dataset lengths (-10.000 
points) for consistent nonlinear measures conflicts with the need for shorter datasct lengths 
(^5.000 points) to provide adequate resolution for the time history generation of trends. 
Consequemly. a ZO-second analysis window was used, as described above, with a 50% overlap 
for an effective tune history resolution of 10 seconds. The nonlinear measures for each 20- 
second analysis window were associated with the time at the center of the analysis window. Le.. 
every tea seconds. Shorter or longer analysis window lengths can be used m proportion to higha 
or lower data sampling rates, as is obvious to those skiUed in the art 

Many characterization tools exist for chaotic data analysis. A subset ofthe tools that 
were found in previous woric to be good measures for EEG data were used in carrying out the 
present invention. TTiese tools include the following: standard statistical measures (minimum, 
maximum, average, absolute average deviation, standard deviation, skewedncss. kurtosis. time 
per cycle); Kohnogorov entropy and entropy spectrum; mmual infomuttion function; maximum 
likelihood correhaion dimension and correUdon dimension spectrum; surrogate generation and 
nonlinearity tests; and nonlinear digital filters (as discussed herein). 

Etmopy. condation dimension, and mutual information were used in the present 
invcntionasnonlinearmeasuresforseizureanalysis. The first minimum in the mutual 
infomiation defines the time scale for generating the return map for EEG dynamics. The return 
map underiies the correlation dimension (measure of dynamic complexity) and emropy (measure 
of dynamic predictabiUty). We have appUed these same measures successfiiUy in analyzing other 
systems. 
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The statistical measures for the present study are obtained by standard methods (6). The 
maximum and minimum are obtained as maximum and minimum (respectively) over the Xj 
values in a time-serial window of w points. The average 00 is given by: 




Eq.6 



The r-th order moment (m.) of the x-data is: 

TTie absolute average deviation (a) provides a robust indicator of the Xi variability (13) and i 
defined as: 



*1 



Eq.8 



An unbiased estimate of the standard deviation (o) is: 



wm. 



w-l 



Eq.9 



An estimate for the skewedness (s) is: 



m. 



Eq.lO 
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An estimate for the kurtosis (k) is: 



III 



» Eq.ll 



The average cycle time (TJ is important as a characteristic time of the nonIii«ar system: 
Y = window length i n timesteps 

I wiiiit^gr of mean erossin2s \ Eq 12 



Tie mmual information function (MIF) is a nonlinear vendon of the (linear) auto- 
cortehition and crossncorrelation functions, and was developed by Fiaser and Swinney (3) 
Mumal information measures the certaimy with which a measurement can be predicted, given the 
outcome ofanotherretatedmeasurement Exan^les of the later include the same EEG cham>el 
at a different time, and another EEG chamiel at the same (or different) time. The MIF indicates 
the average information (in bits) that can be inferred from one measurement about a second 
measurement, and is a function of the time delay (number of time steps) between the 
measurements. The mutual information function also measures the nonlinear time dependent 
correlation in the same signal. For EEG data, infomiation decay in an individual chamiel 

(univariate MIF) indicates local time scale, as the average time lag (V to 

independent of x(t^ and corresponds to the first Oocd) minm 

(3). Fo'«»J««in.aminimumisdefinedastwosuccessivedecieasesinthesi^ 

followedbytwosuccessiveincteasesinsigmd value. Other definitions were tested and found to 

yieldlessconsisteBresults. The MIF. I(Q.S). and system entropy (H) for two measmements(Q 

and S) are defined by: 

/(GhS) = 1(S.Q) = H(Q)*HC^'H(S,Q) Eq. 13 
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Eq. 14 



i 



Eq. 15 



S denotes the whole system that consists of a set of possible messages (measurements for 
the value of s), s„ s^, . . . , s. wi& associated probabilities P$(s,), ?sihif * - 1 Ps(Sii)- Q denotes a 
second system 4at consists of a set of possible messages (measured values wiA a time dday 
relative to the s. values), qi, qj, . . . , qa with associated probabilities ?^qi\ PQ(q2), . • . . PQ(qJ- 
The function PsQ(Sj, qj) denotes the jomt probability of both states occurring simultaneously. If 
the logarithm is taken to the base two, then H is in units of bits. Frascr and Swinney (3) describe 
the details for evaluating I(Q,S), including a sequence of recursive partitions in (Sj, space to 
achieve ndftq^^ft^*^ accuracy with tailoring to the local data structure. 

The maximum-likelihood correlation dimension (D) is based on the early work by Takcos 
(14) with modifications for noise (13): 



\^*cre M is the number of randomly sampled point pairs, r^ is the normalized maximum-norm 
distance between the (randomly chosen) i-j point pairs as defined in Eq. 1 7 (below), and r, is the 
normalized distance (scale lei«th) associated with noise as measured from the time serial data. 
The distances arc normalized with respect to some nominal scale lengtti (Lo), Le. r^ - Lj/Lo and r, 
= LJU> with Lo as a representative scale length (typically a multiple of the absolute average 
deviation). The choice of scale length is a balance between a small scale for sensitivity to local 
dynamics (typically at 1^ < 5a) and avoidance of excessive noise (typically at Lo i a). The 
distances (Ly) are defined by: 




Eq. 16 
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= max |x,,»-x I 



where m is the average number of points per cycle from Eq. 12 Ci-c. m - TJ. Schouten ct al. 
(40) describe the details for evaluating Eqs 16-17 to measure of the number of degrees of 
freedom in a system (e.g^ the number of coupled fiists)nier differential equations to depict the 
5 dynamics). 

The Kohnogorov entropy (K-entropy or simply entropy) measures the rate of infoimation 
loss per unit time, or (alternatively) the degree of predictability. A positive, finite entropy 
generaUy is considered to be a clear demonstration that the time series and its underlying 
dynamics are chaotic An infinite entropy indicates a stochastic, non-detenninistic (totally 

10 unpredictable) phenomenon. For entropy detennination. one begins with two orbits on a chaotic 
attractor that are initially very dose together. Tlie entropy then is estimated from the average 
divergence time for pairs of initiaUy<lose orbits. More precisely, the entropy is obtained from 
. the average time for two points on an attractor to go from an initial separation (L < Lo), to 
become separated by more than a specific distance (L i Lo). The maximum-likeUhood entropy 

15 (K): 



Eq.l8 




£q.l9 



with bi as the number of timesteps for two points, initially within L<Lo, to diverge to Li Lo. Tbe 
work by Schouten et al. (12). and references therein, provide details of the method. Note that the 
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entropy used here is the order-2 Kohnogorov entropy which hereafter is called simply entropy. 

The entropy (K) and correlation dimension (D) usuaUy arc reported in the limit of zero 
scale length. However, EEG data (and all biomedical data) have substantial noise. Consequently, 
the nonlinear measures, K and D, are reported for a finite scale length (Lo) that is slightly larger 
than the noise. Thus, the values of K and D, that are reported here, do not capture the full 
complexity of brain dynamics, i.e., their values arc smaller than expected for the zero-scale- 
length limit. K and D are interpreted as nonlinear statistical indices of finite-scale dynamic 
structure. 



EXAMPLE I 

Hie method and apparatus of tiiis invention are illustrated by the analysis of two datasets. 
Both A»t^>^ are firom the same padent, \Aiio is a 20-year-old female wiA a lifelong history of 
seizures, beginnii^ at age 4 months. The cause of die seizures is not established, aldmi^ 
nemo-imaging soidies (indudii^ computerized tomography and nu«netic resonance) are notmaL 
The seizures are poorly controlled despte treatment widi various combinations of anti-eiMlepdc 
drugs, which at die time of die recordings were Phenytoin, Phenobarbital, and Felbamaie. The 
seizures are partial complex widi some occasions of secondary generalization. During die 
seizure designated Example 1, die patient was sitting up in bed. doing some neurophysiologic 
tesdng. Her EEG shows an activated pattern. She then showed automatisms widi picking 
movements and staring, followed by vocalizations (several seconds of screams). Hyper-extended 
head and neck posturing followed. Her upper extremities became flexed, and dien she showed 
dome activity, involving abduction/adduction at die shoulders and hips. There was ionic 
posniring and clonic activity of aU extremities. "Ibe convulsive movements were associated widi 
high-amplimde EEG waves, involving spikes, polyspikes, and much artifiict activity. As die 
clinical seizure spontaneously termmated, die subject was unresponsive and made loud snoring 
sounds. Then, die brain wave ampUtudes became quite suppressed. The automatisms were 
associated widi polyspike discharges &om die left frontal region. After seizure termination. 



18 



W097/MS24 



PCr/US97/04421 



spike discharges occurred ftom this same region, foUowed by suppressed background waves. 

•Die results from Example I involved four analyses (e-. f-, g., and h-data) on the three 
epochs of cham»el 13 data (epUeptic seizure and post-seizure, non-seizure, and transition from 
non-seizure to seizure). TTiese three epochs of data were provided and analyzed as three separate, 
non-contiguousten-minute data segments. TTie non-seizure data segment ended several hours 
before the transition data, v«ch in turn ended -.1 0 seconds before the start of 
segment The results are combined as one thirty-minute set of plots for the e-data (Fig. 4.1). f- 
data (Fig. 4.2). g-data (Fig. 4.3). and h-data (Fig. 4.4), with gaps in the analysis to indicate where 
one dataset ends and the next dataset begins. In particular, the data from 10-590s is the non- 
seiztire epoch with large amplitude signals at 520.600s from chewing and drinking. The data 
fiom 610-1 190s is the (transition) period immediately before seizure. n« data from 1210.1790s 
includes the seizure and post-ictal phases. The various measures were obtained for every 
analysis-window, and the resulting values were plotted at the center of the 20.second analysis- 
window. 

Each plot in Figs. 4.1 - 4.4 for Example I dispbys the analysis-window^entered measure 

as a solid line. The dashed line (...) in each figure is the average value of the measure (from 
Eq. 6) over an 1 l-point avciaging.window. plotted at the cemral (sixth) point of the averaging 
window. The dotted line (...) in each figure is the corresponding sample standard deviation over 
this 1 1 -point averaging-window (ftom Eq. 9). also plotted at the central poim of the averaging- 
window. TTie scale length (U) was fixed at ~1 .4 times the absolute average deviation, as 
obtained by averaging over the complete non-seizure e-data, and was used in all the analyses as 
the reference scale length for aU three data epochs. Smaller values for this scale length caused 
nmnerical problems in the determinations of the correlation dimension and the Kolmogorov 
entropy; larger values limited the resolution of the nonlinear measures. 

The clinical seizure in Examplcl occurred from 53 to 95s in the (third) seizure epoch 
(1253 to 1 295s in Figs. 4. 1-4.4). Rhythmic convulsions began at 1295s, and post-ictal feanires 
appeared at 1314s. Table 1 shows features in the nonlinear measures that uniquely indicate the 



19 



wo 97/34524 



PCT/US97/04421 



seizure; starred entries {*) denote no clear indicators. Two measures (the peak in Kolmogorov 
entropy of g-data and the peak in correlation dimension of h-data) show epilepsy onset beginning 
-30 seconds before the clinical sei232re. 
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Table 1 Summary of seiznre iodicaton in Example I 
Specific measure e-data 



Time per cycle (TJ 
(timestq)s percycie) 



Entropy (K) 
(bits/second) 



f-data 



g-data 



b-data 



abn^K decrease * • 

T,<200 

125S-1300S 

•bnipiincreaje abnipl incieases abnipt inct«Be abnipt jBcwse 
K>0.063 K>0.008 K>0.056 K>0.016 
1230.1350s 1230.1240s 1220.1340s 1230-1320s 

1270.1275s 

1280.1300$ 



I*Min.inMIF(M,) 
(timesteps) 



Correlation dimension (D) abnpt 



abnqit decrease * 

M,<85 

1290.1310s 



■impt decrease 

M,<20 

1270.1280s 



•ncrease abnipt increase abrupt increase abnipt increase 
D>5 D>5.6 D>32 D>2 

1270.1305s 1340-I345S 1235.1340s 1220.1325s 



EXAMPLE n 

•me dataset designated Example H is fiom the same patient as Example I. who is a 20. 
ycar.,ldfenulewithaUfelong history of seizures as previously described 
the seaure designated Example 0, the patient is lying in bed. awake with her right upper 
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extremity in a flexed posture. The EEC shows spike-wave discharges in tiie left hemisphere, 
after vAdch the brain waves become very sharp, dommated by high frequency activity and 
artifacts. The patient showed eye deviation to the right and some head tunung to the right, with 
head jerking also to the right, but wiftout any usual posturing of the ri^t upper extremity. The 
eye/head turning is preceded by a high-frequency vocalization. After the clinical seizure 
spontaneously terminated, the EEG shows high amplitude wave slowing and subsequent 
amplitude suppression. The subject remained awake during the seizure but was pooriy 
responsive. 

The results from Example II for channel 13 involved four analyses (c-, f-, g-, and h-data) 
on one 23-fninute dataset that included all three epochs of channel 13 data (epileptic seizure, 
non- 

sdzurc, and transition from non-seizure to seizure). The non-seizure period spanned 10-400s* 
The transition period occurred over 4I0-1200s. The seizure began at 1245s and ended at 1290s. 
* The patient was aphasic at 1300s, with head movements and verbalization at 13 15& 

The scale length (Lq) was fixed at -1 .0 times the absolute average deviation, as obtained 
by averaging over the non-seizure e-data, and was used in all the analyses as the reference scale 
length for all three data epochs. The results were obtained as before, and are plotted for e-data 
(Fig. 5.1), f-data, (Fig. 5.2), g-data (Fig. 5.3), and h-data (Fig. 5.4). The seizure is indicated 
cleariy by several unique features in the nonlinear measures, as shown in Table 2; starred entries 
(♦) denote no clear indicators. None of these measures show the epilepsy onset before the 
clinical seizure. 
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Table 2 Summary of seirure indicators in Example II 



Specific measure 



e-data 



f-data 



g«data 



h-data 



Time per cycle (TJ 




• 


abnxpt increases 


abnipt increases 


(timesteps/cycie) 






T,>24 


T.>62 








1255-I275S 


1250.1275s 








1290-1300$ 


1290- 1300s 


Entropy (K) 


sbnipt increase 


iinnifftf III! ■»■■■ A 
au<up& ulCrCaSC 


abrupt incmse 


abrupt increase 


(bits/second) 


K> 0.089 


K> 0.006 


K>0.1 


K> 0.014 




12SS-1290S 


1270.1330s 


125S-1290S 


12S0.1320S 


l«Min.inMIF(M,) 


* 


sbnipt decrease 


• 


• 


(timesteps) 




M,<80 










12S0-1300S 






Correlation dimension (D) 


sbnipt increase 


abnipt mciease 


abnipt increase 


abnipt increase 




D>4 


D>4.5 


D>3 


D>2J 




12SO-1300S 


1280-13005 


1240-1310S 


I240-1315S 



In Example I and Example n, abrupt increases in both the Kolmogorov entropy and 
coirelation dimension indicate both epileptic seizures, although the seizures are clinicaUy very 
different, as described further herein. Abrupt decreases in (TJ number of timesteps per cycle 
(for f-data) and in (M.) the first minimum in the mumal information function (for f- and h-data) 
also indicate the first seizure (Example 1). An abrupt decrease in M, (f.<lata) and abrupt 



increases 
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in Te (g- and h-data) marie the second seizure (Example II). These indicators of seizure are 
summarized in Tables I and 2. The difference in seizure indication (J^ for the two seizures 
implies that care is needed in catalogmg such features. 

It has been sho^ by the foregoing that certain indicative trends such as sudden increases 
and sudden decreases in time serial sequences of nonlinear measures hereinbefore described are 
seizure indicators. 

By comparing indicative trends in non*iinear measures of the patients brain wave data 
with known seizure indicators, it can be clearly determined by whether the indicative trends 
correspond with known seizure indicators whether a seizure is occuning in the patient 

Tables 1 and 2 and Figures 4 and S cleariy show that epileptic seizures are indicated by a 
ntmber of trends in nonlinear measures of brain-wave data including the following: an abrupt 
decrease in the number oftimesteps per cycle for f-data 101; an abrupt increase in tfie number of 
time^teps per cycle for g-data 103; an timspt increase in the number oftimesteps per cycle for h* 
data 105; an abrupt increase in the entropy of e-data 107; an abrupt increase in the entropy off - 
data 109; anabruptincreaseintheentiopy of g-data 111; an abrupt increase in the entropy of h- 
data 113; an abrupt dmra^ in th? mtnimnm m the munifll information function for f-data 
115; an abrupt decrease in the first twmtmum in the mutual information function for h-data 117; 
an abrupt increase in the correlation dimension of c-data 119; an abrupt increase in the 
correlation dimension of f-data 121; an abrupt increase in the correlation dimension of g-data 
123; an abrupt increase in the correlation dimension of h-data 125; and combinations thereof. 

While there has been shown and described what are presently considered the.prcfenred 
embodiments of the invention, it will be obvious to those skilled in the art that various changes 
and modifications can be made therein without departing &am the scope of the inventions 
defined by the appended claims. 
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What is claimed is: 
0) 



(b) 
(c) 
(d) 



A method for automaticaUy detecting an epUeptic seizure in a patient comprising the 
steps of: 

(a) providing at least one channel ofa patient's raw brain wave data. caUedcKiata^ 
selected fiom the group consisting of electroencephalogram dam and 
magnetoncephalogram data; 

separating the e-data mto artifect data, caUed f^iata. and artifect-fiee data, called 
g-daia, while preventing phase distortions in the data; 
processing g-data through a low-pass filter to produce a low^ass-filtered version 
of g-data, called h-data; 

applying at least one nonlinear measure selected fiom the gtoi^ consisting of time 
steps per cycle, kohnogorov entropy, first minimum in the mutual information 
function, and correlation dimension to at least one type of data selected from the 

group consisting of e^ f^ gHlata. and h^ to provide at least one time 
serial sequence of nonlinear measures, fiom which at least one indicative trend 
selected fiom the group consistmg of abrupt increases and abrupt decreases can be 
detennined; 

(e) comparing at least one indicative trend with at least one known seizure indicator, 
and 

determining fiom said comparison whether an epileptic seizure is occurring in the 
patient 



(0 



(2) The method as described m Claim 1 wherein said at least one time serial sequence of 
nonlinear measures is selected fiom the group consisting of: time per wave cycle fore- 
data, time per wave cycle for f^ time per wave cycle for g-data. time wave per cycle 
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for h-data, Kolmogorov entropy for e-daia, Kolmogorov entropy for f-data, Kolmogorov 
entropy for g-data, Kolmogorov entropy for h*data, first minimmn in the mutual 
information function for e-data, first mitn'mnTp in the mutual information function for f- 
data, first minimum in the mutual information function for g-data, first minimum in the 
mutual information function for h-data, correlation dimension for e-data, correlation 
dimension for f-data, correlation dimension for g-data, correlation dimension for h-data, 
and combinations thereof. 

The method as described in Claim 2 wherein said at least one seizure indicator is selected 
fiom the group consisting of: an abrupt decrease in the number of timesteps per cycle for 
f-data; an abrupt increase in the number oftimesteps per cycle for g-data; anabrupt 
increase in the number oftimesteps per cycle for h-data; an abn^ increase in the 
entropy of e-data; an abnqit increase in the entropy of f-data; an abrupt increase in the 
entropy of g-dala; an abrupt increase in the entropy of h*data; an abrupt decrease in the 
first minimum in the mutual information function for f-data; an abrupt decrease in the 
first fn^«*"^""^ in the mutual information function for h-data; an abrupt increase in the 
correlation dimension of e-data; an abmpt increase in the correlafion dimension of f-data; 
an abrupt increase in the correlation dimension of g-data; an abrupt increase in the 
correlation dimension of h-data; and combinations thereof. 

The method as described in Claim 1 wherein the e-data is separated into f-data and g-data 
by use of a zero-phase filter. 

The method as described in Claim 1 wherein the low-pass filter is a standard low-pass 
filter selected from the grotip consisting of second-order, thirdH)rder, and fourth-order 
low-pass filters at firequencies between about 35 Hz and about 60 Hz. 
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The method as described in Claim 5 wherein the low-pass filter is a standaid fourth-order 
low-pass filter at about SO Hz. 

Apparatus for automatically detecting an epileptic seizure in a padent comprising: 

(a) data provision means for providing at least one channel ofa patient's law brain 
wave data called e-data selected fiom the gro^> consisting of 
electroencephalogram data and magnetoencepbalogram H^ tn; 

(b) sq»«»tion means for separating e-data into artifert data, called f-data, and artifiurt- 
free data, called g-data, while preventing phase distortions in die data, 
communicably connected to die data provision means; 

(c) low-pass filter means for filtering g-data to produce a low-pass filtered version of 
g-data, caUed h-data, communicably connected to the separation means; 

(d) appUcation means for applying at least one nonlinear measure selected from tijc 
group consisting of time steps per cycle, Kolmogorov entropy, first minimum in 
the mutual information fimction, and correlation dimension to at least one type of 
data selected fiom the group consisting of e-data, f-data. g-data, and h-data to 
provide at least one time serial sequence of nonlinear measures, ftom which at 
least one indicative trend selected from the group consisting of abrupt increases 
and abrupt decreases can be determined, communicably connected to the low-pass 
filter means; 

(e) comparison means for comparing at least one indicative trend with known siezure 
indicators, communicably connected to the qjplication means; and 

(f) determination means for determining from the comparison whetixer an epUeptic 
seizure is occurring in tiie patient, communicably connected to the comparison 
means. 

nie apparatus as described in Claim 7 wherein said at least one time serial sequence of 
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nonlinear measures is selected from the group consisting of: time per wave cycle for e- 
data, time per wave cycle for f-data, time per wave cycle for g-data, time wave per cycle 
for h-data, Kolmogorov entropy for e-data, Kolmogorov entropy for f-data, Kolmogorov 
entropy for g-data, Kolmogorov entropy for h-data, first minimum in the mutual 
infomiation function for e-data, first minimum in the mutual infonnadon fimction for f« 
data, first minimum in the mutual information function for g-data, first minimum in the 
mutual information fimction for h-data, correlation dimension for e-data, correlation 
dimension for f-data, correlation dimension for g-data, correlation dimension for h-data, 
and combinations thereof. 

(9) The apparatus as described in Claim 8 ixAierein said at least one seizure indicator is 

selected from fte group consisting of: an afanqit decrease in the number of timesteps per 
cycle for f-data; an abrqit increase in &e number oftimesteps per cycle for g-data; an 
abrupt increase in the mmiberoftimesteps per cycle for h<Klata; an abrupt increase in the 
entropy of e-data; an abrupt increase in the entropy of f-data; an abrupt increase in the 
entropy of g-data; an abrupt increase in tiie entropy of h-data; an abrupt decrease in the 
first Tninimiim in the mtitual information function for f-data; an abrupt decrease in the 
first minimum in the mxmial information function for h-data; an abrupt increase in the 
correlation dimension of e-data; an abrupt increase in the correlation dimension of f-data; 
an abrupt increase in the correlation dimension of g-data; an abrupt increase in the 
correlation dimension of h-data; and combinations thereof. 

(10) The apparatus as described in Claim 7 wherein said separation means comprises a zero- 
phase filter. 

(11) The apparatus as described in Claim 10 wherein said zero-phase filter is embodied in a 
programmed integrated circuit semiconductor chip. 
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(12) Tht apparatus as described in Claim 7 wherein said low-pass filter means comprises a 
standard low -pass filter selected ftom the group consisting of second-order, third-order, 
and fourth-order low-pass filters at frequencies between about 35 Hz and about 60 Hz. ' 



(13) 



The apparatus as described in Claim 1 2 wherein said low-pass filter comprises a standard 
fourth-order low-pass filter at about 50 Hz. 



(14) The apparatus as described in Claim 7 further comprising notification means for 

providing notification that a seizure is occurring in the patient, said notification means 
being communicably connected to said determination means. 
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